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We study fractional Josephson effect in a particle-number conserving system consisting of a quasi- 
one-dimensional superconductor coupled to a nanowire or an edge carrying e/m fractional charge 
excitations with m being an odd integer. We show that, due to the topological ground-state degener¬ 
acy in the system, the periodicity of the supercurrent on magnetic flux through the superconducting 
loop is non-trivial which provides a possibility to detect topological phases of matter by the dc 
supercurrent measurement. Using a microscopic model for the nanowire and quasi-one-dimensional 
superconductor, we derived an effective low-energy theory for the system which takes into account 
effects of quantum phase fluctuations. We discuss the stability of the fractional Josephson effect 
with respect to the quantum phase slips in a mesoscopic superconducting ring with a hnite charging 
energy. 
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I. INTRODUCTION 

Josephson effect, a hallmark of the macroscopic quan¬ 
tum coherence, has played a crucial role in the re¬ 
search and applications of superconductivity, both con¬ 
ceptually and practically since its discovery [1-5]. Re¬ 
cently, it has been recognized that Josephson effect 
can also reveal the topological aspects of superconduc¬ 
tivity, most notably the 47r-periodic ac Josephson ef¬ 
fect [6, 7] in one-dimensional class D topological super¬ 
conductors [6-11], as well as in time-reversal-invariant 
generalizations [12, 13]. The doubling of the periodic¬ 
ity is tied to the existence of non-trivial excitations - 
localized Majorana zero-energy modes (MZMs) at the 
opposite ends of a topological superconductor. The pres¬ 
ence of the zero-modes allows for the coherent charge e 
tunneling processes between two superconductors. This, 
in turn, leads to the doubling of the flux periodicity as 
compared to the conventional Josephson effect involving 
charge 2e (Cooper-pair) tunneling. 

In the last hve years there has been a surge of re¬ 
search interest in topological superconductors [9-11, 14- 
28]. Apart from the fundamental importance, the search 
for Majorana zero-modes and other non-Abelian quasi¬ 
particles is fueled by the prospects of topological quan¬ 
tum information processing [29-34]. Recent theoretical 
breakthrough indicating that Majorana-based topolog¬ 
ical phases can be accessed in heterostructures involv¬ 
ing a semiconductor nanowire coupled to a conventional 
s-wave superconductor ]10, 11] has sparked a signifi¬ 
cant experimental activity on this subject [19-27, 35- 
37]. Apart from Ref. [20], most of the aforementioned 
experiments, however, have been focusing on the zero- 
bias peak anomaly associated with presence of the zero- 
energy modes [38-48]. 

Fractional Josephson effect provides perhaps the sim¬ 
plest setup where the non-Abelian nature of the Majo¬ 
rana zero modes are manifested. The stability of this 
effect under various realistic situations (i.e. energy split¬ 
ting, disorder, multiple bands, quasiparticle poisoning) 


have been extensively studied in literature [49-61]. The 
standard approach to understand the Josephson effect is 
based on BCS mean-field theory, where the U(l) parti¬ 
cle number conservation is spontaneously broken. How¬ 
ever, due to the mesoscopic nature of the experimental 
setups in engineering topological superconductors, quan¬ 
tum phase fluctuations may play an important role. In¬ 
deed, charging energy is at the heart of many meso¬ 
scopic superconducting devices such as superconducting 
qubits. Furthermore, charging energy is important for 
topological quantum computing schemes with Majorana 
zero modes [33, 34, 62-67]. There have been several 
works that take into account quantum phase fluctua¬ 
tions in a phenomenological way [55, 57, 68-71]. How¬ 
ever, a microscopic theory of fractional Josephson effect 
is still lacking. In this work, we start from a microscopic 
model of helical nanowires coupled to a fluctuating one¬ 
dimensional s-wave superconductor introduced in Ref. 
[72]. Although there is no long-range superconducting 
order, there are still Majorana zero modes and a related 
topological degeneracy when two or more nanowires are 
coupled to the same quasi-one-dimensional superconduc¬ 
tor (QSC) [72, 73]. We study the Josephson effect within 
this model and find that fractional Josephson effect sur¬ 
vives quantum fluctuations. However, the splitting of 
the ground state degeneracy as well as the hybridization 
between different topological sectors changes and now be¬ 
comes power-law dependent on the system size. In order 
to show that we consider instanton tunneling events be¬ 
tween different topological sectors. We find that phase 
slip events at weak links caused by the spatial inhomo¬ 
geneities (e.g. impurities) in the quasi-one-dimensional 
superconductor are responsible for the power-law decay 
with length of the ground-state energy splitting. Fi¬ 
nally, we calculate the periodicity of the Josephson cur¬ 
rent on magnetic flux through the ring, see Fig. 1 for 
the proposed setup. It has been previously believed that 
one needs to perform an ac measurement to detect Jtt- 
periodicity of the supercurrent [6, 7] which might be quite 
challenging due to parity switching processes (quasiparti- 
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FIG. 1. Proposed setup to study Josephson effect in particle- 
number conserving system, a) The nanowire is coupled to 
a quasi-one-dimensional superconductor forming a loop with 
the circumference L. The four Majorana zero modes 71 , 2 , 3,4 
are located at x = L — I — I', L — I',0,1 respectively. The dis¬ 
tance between 72 and 73 is I' and the distance between 71 and 
72 is Z, same as the distance between 73 and 74 . A magnetic 
flux / pierces through the loop and induces the supercurrent 
flow in the loop. The periodicity of the supercurrent on flux 
allows one to probe topological properties of the system. In 
the case of parafermion proposals, the Luttinger liquid corre¬ 
sponds to, for example, an edge of a two-dimensional Abelian 
fractional quantum Hall system. 


cle poisoning) [9, 10] as well as Landau-Zener transitions 
to the quasiparticle continuum [-57, 74] which together 
put constraints on the appropriate time window for per¬ 
forming the experiment, see below. In this paper, we 
show that by suitably designing the experimental system, 
one can avoid the aforementioned problems and may be 
able to detect 47 r-periodicity in the dc experiments. 

Our approach based on Luttinger liquid formalism 
allows one to extend our theory to the recently pro¬ 
posed systems hosting parafermionic zero modes [75- 
77] . In this case, we consider an edge of two-dimensional 
Abelian fractional quantum Hall system properly coupled 
to QSC. Our main conclusions discussed above for the 
Majorana case remain also valid for parafermions, and we 
show how to probe topological properties of parafermions 
in the dc measurements, see Fig. 1. 

The paper is organized as follows. We begin with the 
qualitative discussion section II where we explain our re¬ 
sults using a simplified toy model which captures some 
main aspects of the fractional Josephson effect. This 
model, however, does not take into account supercon¬ 
ducting quantum phase slips. In order to treat quan¬ 
tum fluctuations properly, one needs to develop a differ¬ 
ent formalism using Abelian bosonisation technique. For 
pedagogical reasons, we first review Josephson effect in 
a Luttinger liquid coupled to a bulk s-wave (Sec. Ill A) 
and p-wave (Sec. IIIB) superconductors, and discuss the 
mechanisms for the change of a periodicity with magnetic 
flux in these two systems. Next, in Sec. IIIC we review 
the fractional Josephson effect in the parafermion sys¬ 
tems by considering an edge of a two-dimensional Abelian 


fractional quantum Hall system at the filling fraction 
V = 1/m coupled to a bulk superconductor. In Sec. IV, 
we present our main results for Josephson effect in a 
number-conserving setup and discuss the dependence of 
the ground-state energy of the system on magnetic flux. 
Finally, the effect of quantum phase fluctuations on the 
flux periodicity of the Josephson current is discussed in 
Sec. V. Some technical details are presented in the Ap¬ 
pendix A. 

II. QUALITATIVE DISCUSSION OF MAIN 
RESULTS 

First, we discuss fractional Josephson effect using the 
simple model involving bulk superconductor with the 
long-range order. Let us consider the setup shown in 
Fig. 1. There are four Majorana zero modes 71 , 2 , 3,4 re¬ 
siding at the ends of the ID topological superconductors, 
and the corresponding effective low-energy Hamiltonian 
reads 

$ 

iJ = zi?,77273 cos —-|-i( 5 ifi 27 i 72 +iJA 347374 -|-iJifi 47 i 74 . 

( 1 ) 

Here Ej is the Jtt Josephson coupling at the junction due 
to the hybridization between Majorana modes 72 and 73 , 
SEij are the energy splittings between Majorana modes 
7 i and 77 which are exponentially small with the distance 
between them SEij oc exp(—Lij/^ij) with being 
the superconducting coherence length in the correspond¬ 
ing segment. It is important to notice that SEij, in par¬ 
ticular 5 Ei 4 , are independent of the magnetic flux due to 
the large superfluid stiffness of the bulk superconductor. 

We now review the physics of the ac Josephson ef¬ 
fect [ 6 , 7]. For pedagogical reasons, it is useful to assume 
that 6 Ei 4 —)■ 0 and 5 Ei 2 = 6 E 34 . In the thermody¬ 
namic limit SE 12 —t 0, the spectrum of Andreev levels is 
given by EEj cos(<i>/ 2 ) where different states correspond 
to even/odd fermion parity of the modes at the junction. 
If the fermion parity is preserved over the time evolution 
of the superconducting phase <i> from 0 to 27r, one can 
see that the system ends up in the excited states, and 
thus the Josephson current through the junction will be 
Jtt periodic. If we take into account finite-size effects, 
one can show that Andreev levels do not quite cross at 
$ = TT due to the exponentially small splitting energy 
SE 12 , see Fig. 2. Nevertheless, provided the evolution of 
the phase <i>(t) = at is fast enough (i.e. a SEx-^, the 
fermion parity of Andreev levels will be approximately 
preserved, and one may still hope to detect Jtt Joseph¬ 
son effect in ac measurements. This phenomenon was 
dubbed as Fractional ac Josephson effect. 

It was shown recently [9, 10, 57, 74] that in order to 
understand the dynamics of the Andreev bound states, 
which is important for the experimental detection of the 
fractional ac Josephson effect, one also needs take into 
account relaxation processes due to quasiparticle poison¬ 
ing as well as Landau-Zener transitions to the quasipar- 
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FIG. 2. Fractional ac Josephson effect. The spectrum of 
the two Andreev bound states corresponding to different 
fermion parity is hybridized at $ = tt by the exponentially- 
small energy splitting SEi 2 . Diabatic passage through the 
avoided level-crossing allows one to follow fixed fermion par¬ 
ity state and to detect the 47 r Josephson effect in ac mea¬ 
surements. The rate of the superconducting phase advance 
is constrained by the parity relaxation (quasiparticle poison¬ 
ing) and Landau-Zener transition rates defined by l/rqp and 
I/tlz, respectively. 


tide continuum, see Fig. 2. The former corresponds to 
parity relaxation processes due to the presence of stray 
non-equilibrium quasiparticles whereas the latter repre¬ 
sents the Landau-Zener transitions during the diabatic 
passage. Roughly speaking, the rate at which the phase 
$ is ramped should be fast compared to the splitting 
energy at the avoided crossing at $ = tt, and fermion 
parity relaxation rate due to non-equilibrium quasiparti¬ 
cles. On the other hand, the sweep rate of cannot be 
too fast, otherwise unwanted Landau-Zener transitions 
into continuum at, for example, $ = 27r would become 
significant. Another relevant timescale th is associated 
with the superconducting phase relaxation dynamics due 
to the resistive environment in which the Josephson junc¬ 
tion is embedded and, as such, depends on the specific 
experimental setup [57, 74]. Overall, experimental ob¬ 
servation of the fractional ac Josephson effect is quite 
challenging because one needs to know a priori the afore¬ 
mentioned timescales. 

We now discuss a new experiment for the detection 
of the fractional Josephson effect using the ground-state 
properties of the system. Let us consider a mesoscopic 
ring shown in Fig. 1, in which the global fermion parity 
p = 71727374 is fixed either by an appreciable charg¬ 
ing energy or by galvanically isolating the nanowire- 
superconductor system so that external electron tunnel¬ 
ing (quasiparticle poisoning) is prohibited energetically. 
Next we assume that the distances L 12 and L 34 are much 
larger than the superconducting coherence length ^ so 
that 5 Ei 2 = SE 34 — 7 > 0. (Otherwise, the flux periodicity 
of the current in the ring is Jtt due to the coherent sin¬ 
gle electron tunneling, similar to the persistent currents 
in non-superconducting mesoscopic rings [78], regardless 


whether the system is in the topologically trivial or non¬ 
trivial phase.) Using the constraint on global fermion 
parity 7174 = —P'f 2 l 3 , the effective two-level Hamilto¬ 
nian can be written in terms of fermion parity at the 
junction 17273 = 1 — 2 n with n = 0 , 1 : 

i7p=i = - [Ej cos I + SEu) ( 2 n - 1 ). ( 2 ) 

Without 6 Ei 4 , the ground-state energy of the system is 
27r-periodic with $. However, when SE 14 is finite, the 
ground state energies at $ = 0 and $ = 27r are clearly 
different, see Fig. 3, and, thus, the supercurrent in the 
ring is not 27r-periodic anymore as discussed in details 
in Sec.V. In other words, we predict that provided the 
charging energy of the ring and the splitting energy SE 14 
are sufficiently large, one should be able to detect Jtt- 
periodic Josephson effect in the dc measurements. We 
note that a crucial requirement of our proposal is that 
global fermion parity is fixed. In the opposite case (i.e. 
global fermion parity is not conserved), the ground state 
becomes 27r-periodic, see Fig. 3c. 



(a) (b) 



(c) 


FIG. 3. Schematic dependence of the many-body energy spec¬ 
trum on magnetic flux (in units of SG flux quantum) for fixed 
overall fermion parity (a) P = 1 and (b) P = — 1, and (c) fluc¬ 
tuating fermion parity. Red (blue) lines correspond to even 
(odd) fermion parity at the junction, (c) The case when 
fermion parity is not conserved. One can notice that the 
ground state energy (solid black line) is 27r-periodic. 

So far we have completely neglected the superconduct¬ 
ing phase slips which are important in mesoscopic struc¬ 
tures. In order to take into account quantum phase fluc¬ 
tuations properly, we will use the bosonisation technique 
which is well-suited for the problem at hand. We will 
show that quantum phase slips affect the level crossing 
between different fermion parity states (i.e. the splitting 
energy E 12 becomes power-law dependent with the sys¬ 
tem size) and, as such, are detrimental for the fractional 
ac Josephson effect. On the other hand, quantum phase 
slips are relatively innocuous in the context of the frac¬ 
tional dc Josephson effect, and our microscopic results 
using the bosonisation technique are in qualitative agree- 
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merit with those obtained using the simple toy model 
above. 


III. JOSEPHSON EFFECT IN LUTTINGER 
LIQUIDS 

In this section we first review Josephson effect in a 
spinful Luttinger liquid coupled to a conventional s-wave 
superconductor [79-82] and obtain the spectrum of An¬ 
dreev states as a function of superconducting phase dif¬ 
ference across the junction. Next, we discuss a spin¬ 
less Luttinger liquid coupled to spinless p-wave super¬ 
conductors [83]. The latter supports Majorana zero- 
energy modes which are ultimately responsible for the 
change of fundamental periodicity of the Josephson cur¬ 
rent. Finally, we will review Josephson effect in more 
exotic structures involving parafermions. 


A. Josephson effect in a Luttinger liquid coupled 
to an s-wave superconductor. 


Let us consider Luttinger liquid coupled to two bulk 
s-wave superconductors with the superconducting phase 
difference $. After integrating out the superconducting 
degrees of freedom, the effective Hamiltonian for the sys¬ 
tem becomes (H = c = 1) 


H = TLnw + J dx {A{x)tjj^{x)tpi{x) + H.c.) (3) 


where m* and /r are the effective mass and chemical po¬ 
tential. One can introduce magnetic field into the Hamil- 
tonian(3) via minimal coupling —idx —t —idx + eA and 
then gauge it away so that it only appears in the super¬ 
conducting pairing potential. Thus, the induced super¬ 
conducting pair potential A(a;) is given by 


A(a;) 


Aq 

cc < 0 

0 

0 < X < 1 

Aoe-*^ 

X > 1 


( 4 ) 


We will assume Aq > 0 in the following. In a ring ge¬ 
ometry, see Fig. 1, the phase difference $ can be related 
to the magnetic flux through the loop 4) = 27r///o with 
/ being the flux piercing the ring and /o = |f is the SC 
flux quantum. 

We now perform standard bosonization procedure for 
spinful fermions using the convention [84]: 


V^r,cr 


^ ^-^lir'Pp-Sp)+<T{r<l>„-e„)] 

\/ 2 Tra 


( 5 ) 


where r = ± and cr = ± for right/left-moving fermion 
with t / i spin, and a the lattice cutoff. The effective 


Hamiltonian now reads 


H= fdxY. ^[K^{dx 0 ^f + 

fJs—p,(T 



dx-^^ cos{V 20 p) 

ZTT Q. 


cos V 2 ip^ 



dx-^^ cos(V 20 p 
2 'Ka 


$) cos 


( 6 ) 


Assuming that Aq is large, the proximity-induced terms 
constrain the values of 0 p and to the minima of the 
cosine potential in the corresponding bulk superconduc¬ 
tors 


0 p{x < 0) = 0 and 0 p{x > 1 ) = 
< 0) = 0 and (pa{x > 1 ) = 


Trip-I-4) 

~7r 


72 ’ 


( 7 ) 

( 8 ) 


where Jp = J 2 r,s and Na = J 2 r,s s^r,s with 

being the zero momentum component of the density op¬ 
erator with spin s and chirality r. Thus, the problem has 
been effectively reduced to solving for the modes of Lut¬ 
tinger liquid subject to the above boundary conditions. 
Note that the allowed integer values for the operators Jp 
and Act have to obey certain constraints (superselection 
rules): (—1)A = (—1)^'^ which immediately follows from 
Eq. (6). 

We now expand the bosonic fields in terms of the nor¬ 
mal modes satisfying the above boundary conditions: 


ipp{x) = + y/lTp ^ ^ cos ^(a„p -k a|,p) 

TT J„ -k $ a; 


0 p{x) = 


n>0 

i 


1 . nnx, j. , 

7wSin^-(a]jCT -Ona) 

— -j= — - —inner “T 

n>0 


n I 

1 . ITTIX 


: Sin 


(^np ^np)- 




n>0 

0 a{x) = \/26»7^ -k -^= ^ ^ cos ^(a„CT + al„). 


Here On and are the annihilation and creation opera¬ 
tors for particle-hole excitations, satisfying the canonical 
commutation relation [amp,a7] = dmndpu- The opera¬ 
tors Lp^p'^ and 6 >ct°^ represent the zero modes of the cor¬ 
responding fields satisfying the commutation relations 
[(/3p°\jp] = i and [0^a\N„] = —i. Using the above 
normal-mode expansion, one can find the energy of the 
system: 


^UVCT,Jpj-^ ^ 1 


T 

Jn -\ - 


p—p,<T k >0 ^ 
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where oj^{k) = Vfj_^ and n^(fc) = (aj,a^). One can no¬ 
tice that the partition function for the system factorizes 
into the product of the zero modes Zq and finite-energy 
excitations Zn- Z = ZqZ^ with only Zq being dependant 
on the flux 

Zo= (10) 

^<y 

Clearly, the sector with odd is gapped out which 
constraints the values of Jp to be even. Thus, the flux- 
dependent ground-state energy of the system becomes 

, , TrKgVg ( $ \ ^ , X 

Finally, one can obtain the expression for Josephson cur¬ 
rent using 


/j($) = 2 




( 12 ) 


One can see that the ground-state energy and the Joseph¬ 
son current through the junction are 27r-periodic which 
is consistent with the previous studies of the Josephson 
effect in s-wave superconductors [79]. 


B. Fractional Josephson effect in topological 
p-wave superconductors. 


Next, we consider a case of topological p-wave super¬ 
conductors and study Josephson effect in the presence of 
Majorana zero modes. Realistically a “spinless nanowire” 
can be engineered in spin-orbit-coupled spin-1/2 quan¬ 
tum wires subject to an external Zeeman field [10, 11]. 
The spinless nanowire is then proximity-coupled to a bulk 
s-wave superconductors at a; < 0 and x > 1. The Hamil¬ 
tonian for the “spinless nanowire” reads: 



-p+iaaydx + Vza:^ 


] tpa'ix), 

aa' 


i7p 


J dx [A{x)ip^ipi + h.c.]. 


(13) 


where a is the strength of the spin-orbit Rashba interac¬ 
tion and is the Zeeman splitting, and superconducting 
pairing is defined in Eq.(4). When chemical potential 
p < Vz, the nanowire is effectively spinless. The elec¬ 
tron tunneling between the NW and the SC leads to the 
proximity effect described by the Hamiltonian iJp. The 
superconducting pairing potential Aq is assumed to be 
a static classical field and quantum fluctuations of the 
superconducting phase are neglected. 

Assuming that Zeeman gap is large, one performs stan¬ 
dard bosonization procedure to find the effective Hamil¬ 
tonian to be equivalent to a spinless nanowire coupled to 


spinless p-wave superconductors: 

H = J dx^[K{d,ef + K-\d,pf] (14) 

— if da: cos 20 + f da; cos(20 — 4)) ) . 

The superconducting phase difference across the junction 
is $. We assume the superconducting pairing potential 
is large (i.e. Ap is relevant and flows to strong coupling 
under the renormalization group flow) and gaps out the 
Luttinger liquid in the region a; < 0 and x > 1. In this 
limit, the values of 9 are constraint to the minima of the 
cosine potential imposing the following boundary condi¬ 
tions for the LL in the region 0 < x < 1: 

$ J- 2TrJ 

9(0) = 0 and e{l) = -. (15) 


where J = Ap — Ap with A^ being the zero momentum 
component of the density operator with chirality r. One 
can expand the bosonic fields in terms of normal modes 
satisfying the boundary conditions: 


ip{x) = + Vk cos + ajj) 

\/n L 

n>0 ^ 


. <I> J- 27r j a: i 1 . Trnx, j. . 

oM = + yf A “">■ 

n>0 


( 16 ) 

Here a„ and aj) are annihilation and creation opera¬ 
tors for particle-hole excitations, satisfying the canonical 
commutation relation [ara,ay\ = 5ran\ is the zero 
mode of the field and is conjugate to J: J] = i. 

After substituting (16) into the effective Hamiltonian for 
the junction, one finds 



One can now easily find the flux dependent part of the 
ground state energy 

ttvK ( 4 ? \ ^ 

We remind that different parity of J = Ap — Ap ac¬ 
tually corresponds to different parity of electron num¬ 
ber operator A = Ap J- Ap. If electron number in the 
junction is conserved, J should be restricted to either 

even or odd sectors, i.e. J = 2m H- ^ ■ Thus, 

the ground state energy as well as the current are 47r- 
periodic: £'g(4>) = £'g(4> J- 47r). On the other hand, 
if there are processes allowing to change the fermion 
parity in the junction, the ground-state energy is 27r- 
periodic. In practice, one should define the time scale as¬ 
sociated with such processes Tp. When t':^ Tp (dc limit), 
the Josephson current in spinless superconductors is 27r- 
periodic and, in this sense, is similar to the Josephson 
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FIG. 4. Ground state energy of a spinless Josephson junction 
(black thick line) as a function of the phase difference. Each 
parabolic curve corresponds to distinct values of J(solid lines 
correspond to even J and dashed ones are odd). 


effect in conventional s-wave superconductors. However, 
if measured at t <C Tp (ac limit), the fundamental period¬ 
icity of the Josephson current is Jtt, and, thus, one could 
distinguish between the topological (spinless p-wave) and 
non-topological (spinful s-wave) junctions. This is why 
this phenomenon in the literature was “coined” as frac¬ 
tional ac Josephson effect [6, 7]. 


subject to the boundary conditions: 

6>(0) = 0 and 0 {l) = ^ + ( 22 ) 

2 m 

Following similar analysis as in the previous section, we 
find the normal mode expansion for the fields 9 and ip is 
given by 


, , /n'\ K 1 TTTIX , j. . 

ip{x) = ip''^> + v — X ^ cos --{an + 

V m ^' vn i 

n>0 '' 


$ J- 2ttJ X 
0 {x) = —^-T + 


2m I ^/Km ^ \/n 

^ n>0 


E l . TTTLX , i , 

^sm—(ajj - an). 


The fields 9{x) and ip{x) should satisfy the commuta¬ 
tion relations (19) and have fundamental periodicity of 
27r. The latter imposes a constraint on the values of J 
requiring that J = 2mk + nj. After substituting above 
expressions for 0{x) and ip[x) into Eq.(20), the ground- 
state energy for the system is given by 




. 27rvK 

= min-— 

k ml 



$ 

47rm 


2 


(23) 


C. Fractional Josephson effect in parafermion 
systems. 


where nj G Z 2 m labels different topological sectors. One 
can see that ground state energy is 47rm periodic Eg{^) — 
Ag($ + 47rm). 


We now review fractional Josephson effect in the pres¬ 
ence of parafermionic zero modes. Instead of a spin¬ 
less Luttinger liquid discussed above in the Majorana 
context, we now consider gapless edge modes of a frac¬ 
tional topological insulator, or the counter-propagating 
edge modes along a trench in a ;/ = — Laughlin state. 
The fundamental excitations in this case are fractional¬ 
ized quasiparticles characterized by the fields ip and 9 
satisfying the following commutation relations: 

[ip{x),9{x')] = -e{x-x'). (19) 

m 

with m > 1. The effective Hamiltonian for the edge reads 

Jdx[K{d,9)^ + K-\d,ip)^]. ( 20 ) 

The physical electron operators in this effective theory 
are given by V'fl/i = thus, the 

proximity-induced superconductivity due to spinless p- 
wave superconductor for a; < 0 and x > I can be taken 
into account using the following effective Hamiltonian: 

Ha = — ( [ da; cos 2m9 + [ da; cos(2m0 — 4)) 

27ra \Jx<o Jx>i / 

( 21 ) 

Assuming that Ap is large, we once again consider Lut¬ 
tinger liquid Hamiltonian confined in the region 0 < x < I 


IV. FRACTIONAL AC JOSEPHSON EFFECT IN 
NUMBER-CONSERVING SYSTEMS 

So far we have included superconductivity at the mean 
field level neglecting quantum phase fluctuations, which 
is appropriate when a nanowire is coupled to bulk super¬ 
conductors. In other words, this approximation implies 
that the particle number is not conserved corresponding 
to the grand canonical ensemble. In mesoscopic struc¬ 
tures, however, particle number fluctuations might be 
suppressed by the charging energy, and, thus, it is inter¬ 
esting to investigate Josephson effect in particle number 
conserving systems. Given that particle number and su¬ 
perconducting phase are conjugate variables, one needs 
to take into account strong quantum fluctuations of the 
superconducting phase in particle number conserving sys¬ 
tems (canonical ensemble). This fact is particularly im¬ 
portant in one-dimensional systems, where U(l) symme¬ 
try can not be spontaneously broken due to the Mermin- 
Wagner theorem. Therefore, in the following we con¬ 
sider the fractional Josephson effect in a model where a 
nanowire is coupled to the quasi-long-range ordered su¬ 
perconductor (QSC) with strongly fluctuating SC phase. 

We now consider an attractive Hubbard model, where 
the spin backscattering caused by the electron interac¬ 
tion is marginally relevant and flows to strong coupling, 
resulting in the formation of the Luther-Emery phase [85] 
with a finite spin gap but no charge gap, and therefore 
can be thought as a one-dimensional analogy of an s-wave 
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superconductor. After the bosonization, the Hamiltonian 
for the superconducting wire reads 

Hsc = H^sc + ( 24 ) 

^sc = (25) 

^sc = [KAd.Oa)^ + K-\d,cl}^f] (26) 

“ [ <ixcos{2V2ip^) 

(27ro)^ 7 

where vp, a and U are the Fermi velocity, the effective 
cutoff length and the interparticle interaction potential, 
respectively. The physics of this quasi-long-range super¬ 
conducting wire can be understood in terms of fluctu¬ 
ating Cooper pairs having algebraically-decaying corre¬ 
lations. The cosine term in Eq. (26) is marginally rele¬ 
vant so that the spin field ipu is pinned to the classical 
minima, opening a spin gap Ag oc e~TET in the QSC. As 
shown below, spin gap is crucial for our construction since 
it prohibits single-electron tunneling to the QSC at low 
energies. Thus, Cooper pair tunneling is the dominant 
tunneling process between the nanowire and the super¬ 
conductor. One can relate the parameters of the above 
model for the Luther-Emery phase to that of a quasi-ID 
superconductor: Kp = 2ny'A^psK and Vp = sjA^p^j k 
with Au,, ps and k being the cross-sectional area of the su¬ 
perconductor, the superfluid stiffness and the compress¬ 
ibility, respectively. 

Our theoretical model also involves the Hamiltonian 
for the nanowire 77 nw (13), in which we once again as¬ 
sume that the Zeeman field is large so that the nanowire 
is in the “spinless regime” with only the lowest band oc¬ 
cupied. The nanowire and the QSC are coupled by the 
single particle tunneling term: 

Ht = t'^ f dx{ijj'lr]a + 'b.-c.). (27) 

a -I 

Here tpa and rja- are electron annihilation operators in the 
nanowire and QSC, respectively. After the bosonization, 
we arrive at following Hamiltonian: 

H = p) + H^sc (^P. T^p) + ^sc (^-: 7^-) + Hr- 

(28) 

Given that single-electron tunneling into the supercon¬ 
ducting wire is suppressed due to the presence of the spin 
gap As, the dominant contribution to the low-energy ef¬ 
fective action comes from pair tunneling. The perturba¬ 
tive expansion in t to second order leads to the following 
imaginary-time action 

6'ph = j dxdTdx'dr' (29) 

a 

■01 (a:, t)iI)^_^{x', T')r]^{x, T)r]-a{x', T')+h.c. . 


Given that the spin held (fa- orders as a result of the last 
term in Eq. 26, the dual field 0a- is strongly disordered, 
and its correlation function decays exponentially: 




\/x^ + VaT^ 


(30) 


This allows us to simplify the action (29) and make a 
local approximation 


5, 


PT' 


(27ra)‘ 


J dr J dx cos ^20 — V 20 p'^ , (31) 


which is valid in the long-time limit \t — t'\ ^ Ag 
Here the Cooper pair tunneling amplitude Ap is given 

The derivation above can be 

apF)‘^ + Vz 

straightforwardly generalized to the case of fractional¬ 
ized Luttinger liquid discussed in Sec. HI C where the 
index m corresponds to a specific edge theory, i.e. m = 1 
represents Majorana case whereas m > 1 corresponds to 
a specific parafermion model. After some algebra, one 
finds that the effective action for the tunneling between 
QSC and LL now reads 



S'pT« — 


A, 


(27ra)‘ 


J dr J dx cos ( 2 m 0 — V^dp'j . (32) 


We also notice that for to > 1, it is not physical to 
consider a finite fractionalized liquid since it exists on 
the edge of a 2D system which has no boundaries. We 
therefore have to induce a distinct gap on the edge to ter¬ 
minate the paired topological regions e.g by a backscat- 
tering term V'pV'L -I- h.c, which becomes cos2m(p after 
bosonization. We assume this is the case for to > 1 in 
the following discussion. 

Finally, one arrives at the effective Hamiltonian for the 
system of interest: 


H = FInw + i^sc + (33) 

i^NW = ^J dx[K{dJ - Af + K-\d^pf] 

^sc = ^ J dx[KpidJp - V 2 Ay + K^^d^PpY] 

iJpT =- — / dxcos(2TO0 — V20p). 

27ra J 

Here we also introduced the vector potential A due to 
the out-of-plane magnetic field, see Fig.l. Henceforth, 
we assume that Cooper-pair tunneling term is large so 
that it “locks” the phase difference between the edge 
modes/nanowire and the QSC 2m0 — ^/20p. [86]. 

Following Ref. 72, we now review the topological de¬ 
generacy in the wires/edge proximity-coupled to QSC. 
Given that the total number of electrons is conserved, 
the minimal setup with topologically protected ground- 
state degeneracy involves four domain walls (i.e. two sep¬ 
arated nanowires) coupled to the same QSC, see Fig.l. 
In the topological regions, 0 = 0-are pinned to 














the classical minima and we denote its value in the first 
and second wire by 0i and 02, respectively. Naively, the 
pair tunneling leads to superficially (2TO)^-fold degener¬ 
ate ground state manifold: ( 0 i, 02 ) = ^(''^■ 1 ,''^■ 2 ),’^■ 1,2 € 
{0,1,..., 2m — 1}. However, one needs to be more careful 
and study the moduli space of the phase variables. In¬ 
deed, since the two wires are coupled to the same QSC, 
we are allowed to make a global gauge transformation 
dp ^ Op + 27r, which leaves the ground state invariant. 
This shows that we need to make the following identifi¬ 
cation: (ni, 712 ) ~ (n-i + k,n 2 + k),k S Z. Therefore, we 
can fix 0i = 0, and different ground states are labeled 
by 02 and we have 2m-fold ground state degeneracy. 

We now study Josephson current for the setup shown 
in Fig. 1. We consider a ring of length L with 0 < x < L, 
and assume that the nanowire covers [0,l]U[L — I — I', L], 
and the QSC covers [0, L — I'] with I' being the length 
of the junction. Using this coordinate system, the four 
Majorana zero modes 71 , 72 , 73,74 in Fig. 1 are located at 
X = L — l — l',L — l',0,l respectively. The vector potential 
can be chosen to be H ^ where / is the magnetic flux 
threading the loop. We first perform the following gauge 
transformation to eliminate the vector potential from the 
Hamiltonian: 


ri<y{x)e 


ip{x) 


'n<y{x) 
iia{x)e~'~9^ 
i>a{x)e 27 - 


L-l-l'<x<L (34) 
0 < a: < ^ 


where $ = 27r///g with / being the magnetic flux pierc¬ 
ing the ring and /o = ff ■ It can be readily checked that 
this transformation is continuous at a; = 0 and, thus, 
fermion operators are single-valued everywhere (i.e. this 
transformation does not introduce a branch cut in the 
fermionic fields). 

However, the pair tunneling terms are affected in a 
nontrivial way and change under gauge transformation 
as 

<») 

Upon standard bosonization, the modified pair tunneling 
terms are given by 

A 

HpT = — —/ da; cos(2m0 — V^Op — $) 

27ra Jq 

A oL — l' 

— —/ da; cos(2m0 — ^/29p) (36) 

27ra 


Assuming that Ap term is large, one can approximate 
the phase difference 9 — by its values at the minima 
of the cosine potential: 




0 L-l-l'<x<L-V 


V2m I X 0 < X <1, 


( 37 ) 


where y = and J G Z. As discussed in Sec. HIC, 

J mod 2m is a conserved quantity, corresponding to dif¬ 
ferent superselection sectors. 

We now calculate the ground state energy and its de¬ 
pendence on the magnetic flux /. It is convenient to inte¬ 
grate out if and cpp fields and write the partition function 
in terms of the imaginary-time effective action. 

pL+l 


S= dr 
Jo 


dx^^^[{dr9f + v‘^{dx9f] 
^JL-i-v 27rz; 


pL-l' 




■vlidjp) 


A 

— —JJ / dx cos(2m0 — V29p — 4)) 

Jq 


(38) 


oL-l' 


- — I dx cos(2m0 — V29p) 

27ra Jl-i-i' 

It is easy to see that the system is gapless in particle- 
number conserving setup contrary to the previous case 
discussed in Sec. HI. Indeed, the combination of the fields 
2m9 + V29p is free to fluctuate. Given that the combina¬ 
tion of the fields 2m9 — \f29p is pinned in the topological 
regions, one can integrate them out. As a result, we can 
impose the following constraint: 

-^drlxOp = dr/x9. (39) 

V2m 

Thus, the problem reduces to that of an inhomogeneous 
single-component Luttinger liquid: 


S = 


[ da;dr-S^. 

/o 27rw(a:) L 


{drO)"^ + v^{x){dx9f 


where the phase field 9 is defined as 

0 < X < L — I' 
L — V < X < L 


9 = 



(40) 


(41) 


The Luttinger parameters K{x) and velocity v{x) are 
given by: 

v+ 0<x<IoyL — I — I' <x<L — I' 
v(x) = L — l — l'<x<L , 

, l<x<L — I — I' 


and 




K(x) = mK 

2 rn?Kn 


with u_|_, IFp being 




/ vK + 2mVpKp 
= y^VpK + 2 mvKp 


(42) 

0<x<l 01 L — I — I'<x<L — I' 

L-l-l'<x<L 

l<x<L — I — V 

(43) 

(44) 

(45) 


KK„ 


Am? 


2 m 


Vn V 
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In order to calculate partition function, one needs to 
specify boundary conditions for the field 6 at x = Q and 
X = L — I'. Since QSC terminates at these points, the 
appropriate boundary conditions are 


where U = Xi — Xi-i and /3 is the inverse temperature. 

After minimizing above expression with respect to 
0 {xi) = Oi,l < i < 3, the minimum of the action is 
given by 


dxOpix) = 0 for X = 0, L — I'. (46) 

The field 9 should be continuous at x = 0, L — Z' so that 
there are no singularities in the effective action. 

9{x~) = 9{x~^) for X = 0,L — I'. (47) 

In terms of the field 9, these boundary conditions trans¬ 
late to 

0(L-/ + O+,r)=0(L-Z-O+,T) 

0(O+,t) = 0(0", r) - X, 

where we used Eq.(7). Thus, the modified field 9 satis¬ 
fies twisted boundary conditions at x = 0. We should 
also include winding numbers of 9 in the imaginary-time 
direction: 


0(x,/3) = 0(x, 0) -I- 27rmM. (49) 


However, one needs to be careful about the values of M. 
Because of the pair tunneling term, M must be an integer 
to preserve the periodicity of 9p. It is easy to see that to 
minimize the action, we can write 


0(x,t) 


27rmM 

/? 


T -I- 0(x), 


(50) 


5™,. = + 


27r/3 
m?L 


loop 


J + 


27r J 


(54) 


where J = 2mk + nj with fc € Z and nj being a num¬ 
ber corresponding to different topological sectors, i.e. for 
m = 1 this number simply denotes fermion parity. The 
total inductance of the loop Tioop = X]i=i i® simply 
given by the sum of all individual inductances (induc¬ 
tors in series rule), and the total capacitance Cioop = 
X]i=i The partition function is obtained by sum¬ 

ming over k: 


Z = ^ exp 

k,M 


27rm^Cioop ^^2 


27r/3 ^ 


(55) 


Now we use the Poisson summation formula to rewrite 
the sum over M: 


5Iexp(- 

M 


27rm^Cioop 

/3 


M' 


= ^exp(- 
N 


7r/5 

2?7i^C7[oop 



(56) 

Therefore, we obtain the ground-state energy for the sys¬ 
tem in a fixed N, J sector (here we restored the units): 


and the action evaluates to 


S = 


27rr7T,^ 


KA,, 




M^ + p [‘'dAAAAAfa 

Jo 


(51) 

Now one can minimize the action (51) for a fixed y. 
Since K(x),v(x) are piecewise constant, the action is 
minimized by piecewise-linear 9. To calculate partition 
function, one needs to find the stationary solution (i.e. 
dr9 = 0) satisfying the following constraints. To sim¬ 
plify the notations, we define xq = 0,xi = Z,X 2 = 
L — I — I',X 3 = L — l',X 4 = L as the locations of the 
interfaces where v and K are discontinuous and the field 
9{xi) are given by 


0(xo) = 0 = 00 

0(xi) = 0i, 1 < i < 3 (52) 

0(X4) = X = 0Ar. 




2 /, 


——^- N + min 

2to"^Cioop Jez nm^Lioop 




2 


(57) 

Eq. (57) is one of the main results of this paper showing 
that for m = 1 the ground state energy is still 47r peri¬ 
odic provided fermion parity is conserved. We can now 
recover the previous results obtained for a bulk supercon¬ 
ductor. Indeed, with the increase of the number of trans¬ 
verse channels A(,h in tbe QSC, quantum fluctuations are 
suppressed. In the limit Kp oc Nch —>■ c», the capaci¬ 
tance and inductance for the loop become Cioop —>■ oo 
and Tioop ~ thus, we recover the results of 

Sec. IIIC. Since J and N commute, the first term is just 
an additive constant corresponding to different particle 
number in the ring which is assumed to be fixed hence¬ 
forth. 


V. TOPOLOGICAL DEGENERACY SPLITTING 


Since we restrict ourselves to the space of piecewise linear 
functions, 0 is then uniquely specified once all the values 
at {xi}?^Q are determined. Thus, the action we want to 
minimize is given by 




N 


i: ^( 9 . 


27r 


k 


(53) 


An important aspect of a topological phase is its 
ground-state degeneracy. In a finite-size system, the de¬ 
generacy is lifted by certain instanton events which con¬ 
nect different topological sectors. In the case of a bulk 
superconductor, such events can be understood in term 
of the overlap of the MZM wave functions. In the inter¬ 
acting system considered here the splitting calculation 
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FIG. 5. Illustration of degeneracy splitting processes, (a) 
electron tunneling between two nanowires through the QSC. 
(b) quasiparticle tunneling across the topological region, (c) 
Electron backscattering in the QSC. 



becomes more subtle. One can identify the processes 
which are present in the BdG formalism and the ones 
that are not, i.e. generated by quantum fluctuations. 
Possible sources of the topological degeneracy splitting 
can be classihed into three types, see Fig. 5: (a) electron 
tunneling between the nanowires through the QSC, (b) 
quasiparticle tunneling across the topological region, (c) 
splitting caused by electron backscattering in the QSC. 

The splitting process (a) is due to the virtual tunnel¬ 
ing of a single electron through the QSC between the 
topological wires (depicted in Fig. 5 as the process (a)). 
We obtain this term in the perturbative expansion of the 
action in the single-fermion tunneling amplitude t: 


Hs 2 


9/2 

^ -A^\x-x'\/vf 

A, 


X cos 



X cos 


vp{x) 


Opix') 

V 2 


+ m0{x) 
— mtp{x') 


(58) 

0 p{.x) _ $■ 

V2 2. 

^p{x') - 


where x, x' are the ends of the nanowire at x = I and 
x' = L — I — I', see Fig.l. Using Eq. (37), one hnds that 
Hs 2 is given by Hs 2 = SEa cos(7r J) with 


2t^ f As a;-a;' 1 

SEa= -y^exp-— 

As V vp 4a: 


■log 


x — x 


■ (59) 


Here we used the fact that ip{x) and ip{x') are pinned 
by the boundaries. For nanowires this is imposed by the 
boundary conditions at the ends of the nanowire whereas 
for parafermion setup the induced backscattering terms 
cos2m(/j effectively pin fields (^(x) and (p{x'), see Refs. 
75-77 for more details. 

The splitting of the ground state degeneracy due to 
the tunneling of a fundamental topological charge can 
be obtained by the instanton calculation corresponding 
to a process that tunnels between two degenerate vacua, 
for example, 1^* “ ^ = 0) and |6» - ^ = ^). This 
process can be intuitively understood as the process of 
a || vortex encircling just the nanowire region, shown 
as process (b) in Fig. 5. The instanton corresponding 
to such an event is homogeneous in space resulting in 
the splitting SEi ^ exp where I is the length 

of the nanowire segment as shown in Fig. 5 and ^ = 


FIG. 6. Ground-state energy dependance on magnetic flux $ 
through the loop for m = 1. Instanton tunneling processes 
SEf, and SEc lead to the hybridization between different topo¬ 
logical sectors whereas SEa process results in even-odd effect. 


v/Ap. Such a process essentially shifts T by 1, and can be 
incorporated into the low energy Hamiltonian as Hsi = 
SEij\J){J -|- 1| -I- h.c.. 

We now discuss electron backscattering effects which 
generate process (c). Consider, for example, electron 
backscattering in the QSC given by the Hamiltonian 

^imp dxS{x — Xi) cos V2ipp{x). (60) 

i 

One can see that ^/2ipp creates a kink of tt in the dual field 
and therefore can be thought as a phase slip event. 

Indeed, a 27r phase slip created in the phase held y/29p 
can be intuitively understood as an vortex tunneling 
across the QSC. Such vortex measures the fermion parity 
in the topological nanowire plus an underlying QSC [87], 
causing a splitting of the ground state degeneracy. Since 
vortex actually measures the charge of the encircled re¬ 
gion, see Fig. 5, such a process is associated with the 
effective charging energy of the enclosed region involv¬ 
ing both the nanowire as well as the QSC. We note that 
impurity scattering in the pair tunneling region is sup¬ 
pressed. Therefore, we consider the effect of impurities 
outside of this region. The details of the splitting calcula¬ 
tion are presented in Appendix A. Here we simply discuss 
our main results. To simplify the instanton calculation, 
we consider the system shown in Fig. 5 with two impuri¬ 
ties at the positions xi and X 2 - Impurity backscattering 
in the QSC leads to the following splitting energy 

He = <5Ac|J)(J-b 1|-bh.c. (61) 

where the energy SEc scales as a power-law of the system 
size. This is an important consequence of quantum fluc¬ 
tuations: topological ground-state degeneracy does not 
scale exponentially as in the case of a bulk superconduc¬ 
tor but rather as a power-law. In the limit Kp ^ K, 
the splitting energy becomes SEc oc viV 2 \xi — X 2 \^~^'’- 
If we replace an impurity at X 2 ^ for example, by the 
hard-wall boundary, the splitting energy is given by 
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FIG. 7. Schematic plot of the supercurrent and inductance 
of the ring as a function of the magnetic flux >1>. Observables 
such as supercurrent and inductance develop Itt periodicity 
once 5Ea is finite. Red (dashed) and blue(solid) lines corre¬ 
spond to SEa = 0 and SEajEj = 0.1, respectively. 

SEc oc vi\xi — which is simply determined by 

the scaling dimension of the cos \/2kpp operator, see Ref. 
72. 

Next, we consider impurity scattering in the nanowire 
described by the following Hamiltonian: 

J dx S{x — Xi) cos 2 irup{x) (62) 

i 

One can notice that an operator acting on the 

ground state, characterized by the field = 9—9pj^/2m, 
shifts 9- by 2 TTm and, therefore, does not induce any 
transitions between degenerate ground states. 

In summary, we have derived low energy theory for 
the superconducting loop structure shown in Fig. 1. The 
corresponding Hamiltonian can be written as 

HeS= ('^+^) +dEa{-iy^ \J){J\ 

+ {SEi) + SEc)\J){J + 1| + h.c., (63) 

2 f ^ 

where J = 2mk + nj and Ej = —. The spectrum 

of the system is shown in Fig. 6. One can see that there 
are two types of non-commuting splitting terms. The 
processes described by SEf, and 5Ec splitting energies re¬ 
sult in the hybridization between different J and J + 1 
sectors and, thus, open a gap at the avoided level cross¬ 
ings. This is not surprising since fractional excitations 
have an associated charge e/m, and thus charging en¬ 
ergy inducing vortex tunneling can distinguish between 
different topological sectors and couple them. On the 
other hand, we have a process causing the splitting en¬ 
ergy 5Ea which lifts the degeneracy between even and 
odd J-sectors. In particular, in m = 1 case the split¬ 
ting energy SEa distinguishes between even and odd par¬ 
ity sectors and restores 47r periodicity of the fractional 
Josephson effect, see Fig. 6. The results obtained via 
the bosonization analysis of a microscopic Hamiltonian 
agrees well at the qualitative level with the simple effec¬ 
tive Hamiltonian (1) subject to the total fermion parity 
constraint, see Sec.H. 


VI. CONCLUSIONS 

In this paper we study Josephson effect in mesoscopic 
superconducting structures and focus on the contribution 
to the supercurrent due to the presence of Majorana zero- 
energy modes or more exotic parafermionic modes. In the 
former case, coherent single electron tunneling across the 
Josephson junction is allowed due to ground state degen¬ 
eracy and results in the 47r-periodic (fractional) contribu¬ 
tion to the Josephson current. Previously, it was believed 
that fractional Josephson effect may only be accessed in 
ac measurements [7] which might be quite challenging in 
realistic experimental settings [74]. Here we show that 
one might be able to detect this anomalous periodicity 
with magnetic flux in mesoscopic superconducting rings 
using dc measurements in which the flux through the 
ring is changed adiabatically, and, therefore, to avoid the 
aforementioned issues. 

The system we consider here consists of a quasi- 
one-dimensional superconductor (no long-range order) 
coupled to a nanowire or an edge carrying e/m frac¬ 
tional charge excitations with m being an odd integer. 
Using bosonization technique and instanton analysis, 
we derive an effective low-energy Hamiltonian for the 
system by taking into account instanton events which 
couples different topological sectors, see Eq. (63). 
We show that quantum phase fluctuations lead to a 
power-law dependence of the splitting energy with the 
distance. This is to be contrasted with the splitting 
energy in long-range-ordered superconductors which 
scales exponentially with the system size. We calculate 
the periodicity of the supercurrent on magnetic flux 
piercing the superconducting loop and show that, by 
properly designing mesoscopic ring, one should be able 
to measure 47r-periodic component of the supercurrent 
in dc experiments. We believe that our proposal is 
within the experimental reach and provide here detailed 
recipe for measuring fractional dc Josephson effect in 
mesoscopic systems: 

• fabricate a mesoscopic ring with approximate di¬ 
mensions of lOOnm in diameter and 10/rm circum¬ 
ference such that the charging energy of the ring 
is much larger than the dilution fridge tempera¬ 
ture. The circumference of the ring should be much 
larger than the superconducting coherence length 
which is estimated to be of the order of lOOnm (i.e. 
Ui 2 ,U 34 » 0: see Fig.l. 

• bring the outer ends of a semiconducting nanowire 
within the superconducting coherence length f so 
that there is an appreciable splitting energy SE 14 
(for Li 4 < the splitting energy SE 14 can be of 
the order of p-wave gap). 

• measure the supercurrent in the ring using, for ex¬ 
ample, SQUID [88] or torque [89, 90] magnetom- 
etry technique. We estimated the supercurrent in 
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the ring to be of the order of lOnA which should 
be within experimental reach. 

• the semiconductor nanowire should be designed in 
such a way so that one can drive topological quan¬ 
tum phase transition with an in-plane magnetic 
field (for example, 17-shape nanowire). In the triv¬ 
ial phase (e.g. at zero magnetic field B = 0), 
the ground state energy exhibits 27r-periodic de¬ 
pendence on the external flux whereas in the non¬ 


trivial phase the flux periodicity should be 47r. 
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Appendix A: Degeneracy splitting due to quantum phase slips 

In this Appendix, we calculate the topological degeneracy splitting due to quantum phase slips. To be concrete, 
we consider quantum phase slips generated by impurity backscattering. For simplicity, we consider the geometry 
shown in Fig. 5 where the impurities are located at the positions Xi and X 2 , just outside the pair-hopping region. For 
simplicity, we consider m = 1 case. The corresponding Hamiltonian then reads 


-ffimp = Vl cos V2ipp{xi) + V2 cos V 2 ipp{x 2 )- 


(Al) 


We now compute the amplitude M for the instanton tunneling between the minima of cos(20 — \/2Bp) in the interval 
X\ < X < X2 '■ 


M = ( 0 - %=TT 


^-HT 


0-^=0 

V 2 


(A2) 


T-»-o 


where H is the total Hamiltonian including i7imp- This amplitude can be calculated using a path integral with the 
appropriate boundary conditions: 


Z= J 'D^'DOSi9-{x,O)-Tr{2ki + l))S{0-ix,T)-2TTk2)e-^‘^^'^^^^+^^-^^^ 


(A3) 


ki 

Op{x,T) 


where 9-{x,t) = 9{x,t) - 


and 


C = -d^ifpdrBp + -d^ipdrB + — [K + K{da;9)'^] + ^[K + Kp{da,9pf] 

TT 7T ZTT ZTT ^ 

— ^^Q{x — xi)Q{x 2 — x) cos[20(a;) — V^Bplx)], (A4) 

27ra 

mp = [ui(5(a: — xj") -I- V2S{x — xj)] cosV2(pp{x,T). (A5) 


Here the pairing tunneling term is non-zero in the interval xi < x < X 2 which defines the topological wire segment. 
In the limit Kp > 2 when backscattering is irrelevant, we can calculate partition function perturbatively in Vi 


Z = y VipV'Be _ J DipjyQ ^ J dridr 2 cos \/2(/3p(X;^ , n) cos •\/2(/?p(xJ, r 2 ) H—e 

Here VB denotes integration with the 6 function constraints, see Eq.(A3). Using the identity 


— f drdxCo 


cosV2pp{x^,Ti)cosV2pp{xt,T2) = j 


Si,S2=±l 


one finds that the effective action for ipp is given by 


(A6) 

(A7) 


5[i^p] = / dxdr 


-da;P>pdr9p + 

" ZTT 


-I- iV 2 [si(^p(xi , n) -I- S2(Ppix^, T2)] . 


(AS) 
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After integrating out the effective action becomes 


= 


d^q 

(2^ 

d^q 




wfc<^p(q)6»p(-q) + v^^‘/?p(q)‘/?p(-q) + + S2e*'i'’"")v3p(q) 

zttA, 


^^^'^p(q)'^p(-q) + i f^^p(-q) + + 72526*'’'^" ) </5p(q) 


(A9) 


where q = (fc, ic), x = (a;, r) and the inner product q • x = fcx — lot. We can now integrate out Lpp to find 


5[0p] = - 


d^q irKp 
(27r)2 2vpk^ 


—0p(q) + + V 2 s 2 e-^^-^ 

7T 


K„ 


d2q 
2TTVp J (27r)2 


000 piq) 


V^TTSl 


2q xi 


V^7rS2 


2qX2 


Kp 

2 'KVn 


drdx 


drOp + V2Trsi<d{x — Xi)6{t — ti) + •\/27rs20(a; — a;2)i5(T — T 2 ) 


n 2 


(AlO) 


Next, we integrate out ip and obtain the following effective action: 


Cs,sA 0 p, 0 ] = ^ 


(j0T0p + •\/27rsi0(a:i — x)S{t — ti) + V2Trs2Q{x2 — x)6{t — T2)^ + Vp{dx0p)^ 


K 

2 tt 


-{drO)"^ + v{dx0f 


— ^^Q{x — xi)Q{x 2 — x) cos [20(a:) — V29p{x)\ . (All) 


Combining all the terms together, the partition function now reads 


Z = 


J vOpV'e 




(A12) 


The calculation of the first term / 'D'9'D9pe~^° reproduces the splitting energy 5Ei, see Sec.V. We will focus here on 
the second term and calculate the contribution of the classical field configuration minimizing the action One 

can notice that the the quadratic action of 9p contains S function in r. Therefore, in order to get a finite action 9p 
field must have a discontinuity at ti and T 2 . Indeed, let us write 9p as 

9p = 9p — A{x),A{x) = y/2'K[siQ{xi — x)0(t — n) + S20(x2 — x)0(t — T 2 )], (A13) 


where 9p{x, r) is a now smooth field as far as time dependence is concerned. Notice that in doing so we have introduced 
a jump in the spatial profile of 9p at Xi and X 2 . The discontinuity at these points has to be carefully taken into account 
by considering an inhomogeneous problem since the pairing field also has jumps at a;i and X 2 . However, in the limit 
when the length of the topological region is large, the bulk energy gives dominant contribution and thus the boundary 
effects can be ignored. 

Next, we rewrite the action using the new fields: 


Cs,s^[9p,9] —{dr9p)‘^ + Vp{da;9p-dxA)'^ 

Ztt 


Vn 


K 

27r 


-'p) \ ^p\^x^p 

^-{dr9f + v{dg,9f 

V 


^^0(X — Xi) 0 {x 2 — x) COS \29(x) — . 

27ra 


(A14) 


In the domain xi < x < X 2 , the combination 9 — 9pj'/2 is pinned and one can use the relation d^/rd = dx/r9pj\^ to 
simplify the calculation. Thus, within this space of field configurations, the corresponding partition function can be 
evaluated exactly 


2 'ViV 2 — 


V 1 V 2 


dri dT2 


V9pe 


—Seii 


5eff[0] = 


1 

drdx — 
2 tt 


K. 


K\ 


vK , 


+ VpKp{d.9p - d^AY + —{d,9p) 


(A15) 

(A16) 
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We can now evaluate the path integral and calculate the dependence of this term on L = X 2 — xi. Let us look closer 
at the term {dxd — A)^: 


dxdr [dxOp — dxAf = / dr / dx 


{dxOpf - 2\/2TT{dx0p){xi,T)<d{T-Ti) + 2\/2TT{dx0p){x2,T)Q{T-T2) + {dxAf 


(A17) 

For our purpose, it is sufficient to consider large L limit and neglect the boundary conditions for 9 which does not 
affect the scaling of the action with L. Therefore, we neglect the divergent boundary term {OxA)"^ in Eq.(A17) which 
is simply an artifact of our approximation. The effective partition function in the momentum space becomes 


^eff = 


dfcdw r 1 / Kp K\ 2 

(27r)2 [ 27r \Vp 2vJ 




|dp(q)P + ^/2vpKp^{s^A^-^-+S2e^-^-)~ep{ci)^ . (A18) 


It is clear that the correlation function is vanishing for si = S 2 - Therefore, in the following we set si = 1,S2 = ~1- 
After integrating out 9 field, one finds 


V9. 


g “Seft^ gjj-p I _ 


dfcdw 


viv+K^pk 


21,2 


Att 


Kj^uA(uA + 


fc2) 


|giq.X2 _e*q-xi|2 


= Cl exp — 


P P In 

2 vlK+ 


vX\ti-T 2\^ + \Xi-X2\ 


(A19) 


where a is a UV cutoff and Ci is the numerical prefactor. Finally, the contribution to the partition function reads 


-Z.Si.So = 


V 1 V 2 C 1 


dTidr2 


In - T 2|2 + |a;i - a: 2 p 


ViV 2 Ta 


y+ 


\Xi - X2\ 


--1 


(A20) 


Following standard calculation, see Ref. [91], in order to obtain the energy splitting, we need to take into account 
multiple instanton processes corresponding to multiple insertions of jumps of 9p. In the end, one finds that the energy 
splitting is given by 


SE 


ViV 2 a 



(A21) 


and is power-law dependent on the system size. In the limit Kp ^ K, v+ ~ Vp and K+ « Kp, so the splitting energy 
becomes 6 E oc Note that the boundary of a QSC can be represented as a strong impurity, say, at X 2 which 

effectively cuts off the superconductor at this point. In this case, the fluctuations of the phase ipp at X 2 are suppressed, 
and the splitting is given by scaling dimension of C0S'/2ipp{xi) [72], i.e. SE oc \xi — X 2 \^~^'’^'^ at Kp ^ K. 
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